Anisotropic Coherence properties in a trapped quasi2D dipolar gas 
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We consider a trapped quasi2D dipolar Bose Einstein condensate (q2D DBEC) with a polariza- 
tion tilted into the plane of motion. We show that by tilting the polarization axis, the coherence 
properties are anisotropic. Such a system will have density fluctuations whose amplitude depend on 
their trap location. Additionally, interference contrast will also be anisotropic despite an isotropic 
density profile. The anisotropy is related to a roton like mode that becomes unstable and supports 
local collapse along the polarization axis. 
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Coherence properties underlie collective quantum be- 
havior. An example of a need to further understand co- 
herence is in high T c superconductors pQ , where the origin 
of the phase coherence leading to superconducting be- 
havior is still mysterious. Such solid state systems where 
this occurs are complex, making the study of coherence 
properties a challenge. One avenue to study coherence of 
quantum systems more clearly is in ultracold gases, which 
have offered direct insight and unprecedented control to 
the study of quantum collective behavior. For example, 
two macroscopic Bose Einstein condensates (BECs) were 
interfered displaying phase coherence [2J. More rehned 
experiments have directly detected higher order correla- 
tions [3] which cannot be easily accessed by electronic sys- 
tems. Additionally, reduced dimensional ultracold gases 
have been the focus of intense interest because fluctu- 
ations are enhanced in such geometries. Phase coher- 
ence has been studied at the Berezinskii Kosterlitz Thou- 
less (BKT) transition [H [5] where two independent q2D 
ultracold gases were interfered, revealing thermally ac- 
tivated phase defects are the origin of decoherence [BJ. 
The BKT transition is a transition where the phase co- 
herence changes character, and insights into controlling 
such coherence properties are of general importance. In 
this paper we show the coherence properties of a trapped 
dipolar gas can be controlled. 

The progress of experiments to probe ultracold gases 
is remarkable. Recent advances in imaging q2D ultra- 
cold gases have been used to detect density fluctuations 
[?HS]- Additionally, ultracold atomic systems in one di- 
mensional systems have also amazingly been used to 
study quantum fluctuations directly 10J and other co- 
herence properties [111114) . For these examples of ultra- 
cold gases, the interaction between the constituents are 
short range. There are ultracold dipolar systems which 
have anisotropic, long range interactions, which offer an 
entirely new avenue to study quantum correlations and 
collective behavior. Amazing experimental progress has 
been made with the strongly magnetic atoms: Chromium 
(Cr), Dysprosium (Dy), and Erbium (Er). All three have 
been Bose condensed and displayed strong dipolar ef- 
fects [I5HT7] . The study of dipoles in reduced geome- 
tries have begun with Cr, where it has been put into a 
one dimensional lattice, forming a series of coupled q2D 



dipolar BECs [T5]. Additionally, polar molecules have 
been loaded into q2D geometries to control the rate of 
their chemical reaction |19j . 

Correlations of ultracold dipolar gases have been stud- 
ied with a perpendicular polarization for weakly interact- 
ing gases near the roton instability [2"0H2"2"] and strongly 
interacting gases, e.g. 23J. Using the direction of the 
polarization axis as a means to control the collective be- 
havior and coherence properties, a q2D dipolar gas has 
been studied in the many body dipolar systems at the 
mean field level [2H US] an d for 2D dipolar scattering 
[26] . Additional theories have investigated a tilted polar- 
ization axis in a q2D dipolar fermi gas [2"T] . 

In this work, we go beyond mean field theory and ex- 
amine the correlations of a trapped q2D DBEC at fi- 
nite temperature with a partially tilted polarization into 
the plane of motion. This situation leads to anisotropic 
coherence properties, in which the gas maintains coher- 
ence along the direction of polarization, but more rapidly 
loses coherence perpendicular to the polarization direc- 
tion. This anisotropy is evident in both the phase and 
density-density coherence properties. This emerges be- 
cause of a roton-like mode in the excitation spectrum 
and is strongly anisotropic in character. We consider the 
implications of this theory through interference experi- 
ments and density fluctuations, both of which will ex- 
hibit anisotropic correlated behavior. We also calculate 
the compressibility of the dipolar gas. 

To study the finite temperature coherence properties, 
we employ the Hartree Fock Bogoliubov method within 
the Popov approximation (HFBP) with nonlocal interac- 
tions [351 I2H] • The HFB breaks the wavefunction into a 
condensate and thermal component: ^> = [s/No^o (p) + 
9(p)] where we have replaced ao — > VNo- We use the 
Bogoliubov transformation: 8(p) — '^2 1 \u 1 {p)a 1 e~' lu '~< t — 
vt i {p)a* 1 e lul '' t ] where a 7 (a*) is the bosonic annihilation 
(creation) operator for the j th quasiparticle. The HFBP 
solves a non-local, generalized Gross-Pitaevskii equation 
and a nonlocal Bogoliubov de Gennes set of equations 
from which the eigenvalues and vectors (quasiparticles) 
are self-consistently obtained [2"5] . 

To evaluate the interaction, we assume a dipole mo- 
ment of the form: d — d[x sin(a) + z cos(a)] where a 
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is the angle between z and d, and we assume a single 
transverse Gaussian wavefunction of width l z is occupied 
in the z direction, where ij = ^Jhju^ra and £Jj is the 
trapping frequency. This requires /i,T <C Sw 2 , and this 
leads to a simple form of the interaction. We use an 
isotropic trap (x,y) with l z /l p — 0.15 {uj z /lo p =4,AA) 1 and 
an interaction strength of gd = 0.025 and g — 0, where 
g c i = cP/\/ 1 2/kI z and g — \f%iajl z with a being the 3D 
s-wave scattering length (a <C l z ). The polarization an- 
gle is set to a/n = 0.25, this is beyond the critical angle. 
The critical angle is when the interaction is zero along 
the x axis, so sin(a c ) = l/\/3- For a > a c ~ 0.27T there 
is an attractive region to the interaction. This value of a 
is picked to optimize the collapse with respect to parti- 
cle number. If it were steeper, the critical number would 
be less, and the depletion would not be large enough 
to observe the desired behavior. In addition, if a were 
smaller then roton the mode would require many more 
particles to collapse, and the collapses would be more 
3D in nature. For a/ir = 0.25, the system is well de- 
scribed by the q2D formalism. For these parameters, the 
trapping potential for Dy [Er] is (w p , u) z )/2-k— (12,533) 
[(48,2133)] Hz. In a trapped ideal 2D gas, Bose conden- 
sation can occur [30] , and the the critical temperature is 
at T c /huj = VM/tt ~ 0.78VN. This is for the pure 2D 
trapped gas, we only use it as a reference temperature. 
Throughout this paper we work in trap units. It is impor- 
tant to mention that the examples given here are weakly 
interacting; the HFBP method cannot accurately handle 
strongly interacting dipolar gases (// < 10Hlj p ) 29J. This 
criteria is from the deviation of the Kohn mode, or slosh 
mode, from lhw p . If the interaction is too strong then 
the static thermal cloud alters the slosh mode of the con- 
densate; this is a numerical signature that the HFBP has 
broken down. This restricted the work from entering the 
Thomas Fermi regime where l p > £, where £ is the heal- 
ing length (y^mfi/h 2 ). Because of this, trapping effects 
are always important in this study. 



Once gd, g, l z , T, a and the condensate number, Nq, 
are selected, the HFBP calculation finds: the number of 
thermal atoms, iV; the chemical potential, /Lt; the con- 
densate wavefunction, 0o; and the quasiparticle wave- 
functions, Uy and v-y. From these wavefunctions, we 
construct the correlation functions: the nonlocal ther- 
mal correlation function, h(p, p 1 ) = J2- y [ u j(P) u 7(p f ) + 
v*(p)v~ r (p')]N 1 + v*(p*)v~,(j?) with N 7 being the Bose- 
Einstein occupation, and the condensate correlation 
function, no(p, ff), given by Nq(/)q(x)4>o(x / ). The to- 
tal local density is n(p) = rio(p) + h(p) where e.g. 
n o(p) — n o{p,P)- The HFBP method uses these correla- 
tion functions in the self-consistent calculation, and they 
can be related to the more common g\ and gi correlation 



functions [31] : 

si(p;/5") = 



<**(pW)> n {pj)+n{pj) 



v / n(p)n{p') v / n(p)n(p') 

(vir(flttV)vi/(pW)) 

n(p)n(p') 



(i) 

(2) 



J\2 



1 +9i(p,P ) 



n{p)n(pT)' 



gi details phase correlations, for example it determines 
the fringe contrast if a BEC were split and interfered. 
gi is the density-density correlation function, and is re- 
lated to density fluctuations and the static structure fac- 
tor 

We study a system with a temperature of T/Tc = 0.5. 
At this temperature, the condensate fraction is No/N > 
0.57 for all N considered. This is well within the region 
of validity for the HFBP. This temperature thermally 
populates the quasiparticles, so that the anisotropy in 
the coherence properties is significant. Furthermore, this 
temperature makes the density more isotropic. Values 
of T/Tq between 0.4 and 0.6 produces similar results. 
Raising the temperature would lead to populating more 
thermal modes, and would make correlations isotropic, 
reflecting the nature of the trap (T/Tq ~ 1). For lower 
temperatures, the correlations are determined by the con- 
densate's properties alone, as the depletion and thermal 
population become negligible. A pure condensate has to- 
tal phase coherence, it is the thermal fraction that leads 
to reduced coherence. 

We now look at the correlation functions, g±(p; 0) and 
.92 (p] 0). These quantify the coherence of the system from 
the origin to p. To illustrate the anisotropic coherence of 
the gas, we plot them along the polarization axis and in 
the perpendicular direction in figure [l] 

Fig. [I] (a) shows the total density and Fig. [I] (b) shows 
gi(p; 0) along the x (solid line) and y (dashed) axes. The 
important feature is that g\ quickly decreases in the y 
direction, but only steadily decreases its value in the x 
direction or polarization direction. We have shown the 
system for three different number of total particles 2000 
(blue), 3000 (red), and 4000 (black), from top to bottom 
at the temperature T/Tq = 0.5. Varying the number 
is like varying the interaction strength. The density is 
isotropic near the center of the trap, even when the gi is 
anisotropic. Near the edge of the trap the total density 
becomes anisotropic and is more extended in the x axis. 

A recent Bragg spectroscopy experiment observed an 
anisotropic speed of sound for a 3D dipolar gas [33]. A 
similar experiment for a q2D DBEC could be used to 
probe these predicted anisotropic coherence properties. 
If a DBEC were split into two equal clouds such that 
one is moved by a with a relative phase of S(p), then 
the total density would be n(p) — 0.5(n(p) + n(p + a)) + 
\J n(p)n(p + a)gi (p, p + a) cos[S f (p)]. The fringe contrast 
will depend on the angle between a and d because of the 
anisotropy in g\ (p, p + a). To maximize the difference in 
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FIG. 1: (a) The total density, (b) gi(a;,0;0,0) (solid line) 
and gi (0, y; 0, 0) (dashed) , and (c) gi (x, 0; 0, 0) (solid line) and 
g 2 (0,y;0,0) (dashed) at T/T c =0.50 and for N =2000 (blue), 
3000 (red), and 4000 (black) from top to bottom. The to- 
tal density is isotropic in the center and yet the correlation 
function is significantly anisotropic at short range. 



fringe contrast, the phase difference should be ir/2 when 
a is near 3/4 l p . This signal is only a few percent and 
would be challenging to measure. 

We now move on to the density-density correlation 
function, g 2 (p; 0) . In Fig. [I] (c) , 172 this is plotted along 
the x (solid line) and y (dashed) axis. The results are 
similar to those of g\ , the extent of the coherence is fur- 
ther in the x direction than the y. The anisotropy is more 
significant in this correlation function. The fact that the 
52 goes below one is a trap effect and is most signifi- 
cant when correlating with the trap center [HI [12] . This 
correlation function is important for density fluctuations, 
this will be discussed below. 

To understand the nature of the anisotropy we look 
at the Bogoliubov de Gennes excitation spectrum as a 
function of particle number, this is shown in Fig. [2] (a). 
The spectrum goes soft or an eigenvalue becomes com- 
plex implying the system is unstable when there are more 
than 4850 particles. Additionally, we show contour plots 
of the wavefunctions of the two modes which go soft (it 7 ) 
in Fig. [2](b,c), and the dashed black lines are y/rio- The 
contours at 0.25, 0.5, and 0.75 of the maximum values of 
the individual wavefunctions, and for it 7 the shaded re- 
gions signify a region less than zero. These quasiparticle 
modes start with an energy at or above Ahuj p . Then as 
the particle number is increased, their energy decreases. 
In a homogeneous q2D DBEC, the roton is a local min- 
imum in the excitation spectrum, which occurs near the 
k p l z ~ 1 [31] ■ If the spectrum goes soft, the system lo- 
cally collapses into plane waves which form dense stripes 
along the polarization axis. 

The Bogoliubov de Gennes excitation spectrum has 
been analyzed in detail in Ref. [53] for a dipolar gas 
with a = or a cylindrically symmetric system. In that 



case there are degenerate modes with azimuthal symme- 
try ±m. In the present case this degeneracy has been 
broken by the tilted polarization axis, leading to slightly 
different energies depending on the excitation shape rel- 
ative to the x axis. This leads to band like structures 
forming in Fig. ^ (a). 

The mode in [2^ (b) makes a significant contribution 
to the correlation function for two reasons: first, it is 
a low energy mode, and is more thermally populated. 
Second, the correlation function, gi{p '■ 0), has a factor of 
it 7 (0) and u 7 (0) in it. Only a few low energy modes have 
a significant contribution at the origin. The breathing 
mode just below 2Hui and the roton mode in [2] (b) are 
two low energy modes with a maximum in the center of 
the trap. Other low energy modes have small or zero 
amplitude in the center of the cloud because symmetry. 

In the trapped case, the finite extent of the gas alters 
the nature of the "roton" mode 35 [371 . For a = 0, the 
gas become 3D in nature, p, ~ hu z [22] . When the gas 
becomes unstable, it locally collapses and forms regions 
of high density. A good demonstration of this is given 
in Ref. [38]. For a/ir ~ 0.25, the stability of the sys- 
tem displays an interplay of the interaction and the trap 
geometry, while maintaining p -C hio z . The collapse is 
local and forms a series density regions parallel to the 
polarization. In Fig. [2] (b,c), we see this character in the 
quasiparticles, which have an extended amplitude along 
the polarization axis and a series of nodes in the perpen- 
dicular direction. It is this character of the quasiparticles 
that leads to anisotropic character to the correlations. 
The nodes in the roton modes lead to destructive inter- 
ference and reduce correlation in that direction. 

To further study the implications of the anisotropic 
correlations, we consider the density fluctuations of 
the gas by calculating the compressibility: kTn{p) = 
(Jdff6n(p)Sn(ff)) =Jdffn(p)n(ff)g 2 (p;ff)-n(p)(N-l) 
where 5n(p) = ^f*(p)^(p) — n(p) and this is the density 
fluctuations of the system. The compressibility and num- 
ber fluctuations have been measured in BEC experiments 
[H [5] ■ We present the compressibility for a trapped dipo- 
lar gas in Fig. [| (b) at T/T = 0.5 for 2000 (blue), 3000 
(red), and 4000 (black) in the x (solid) and y (dashed). 
These results show that the density fluctuations roughly 
follow the anisotropic density profile. Another way to 
state this is: that depending on the local value of the den- 
sity, the compressibility is determined almost uniquely. 

In this work, we found that the correlations are 
anisotropic in a q2D DBEC with a polarization axis tilted 
into the plane of motion. We studied the DBEC with the 
HFBP at at temperature of T/T c ~ 0.5. At this temper- 
ature, quasiparticles are thermally populated and total 
density is nearly isotropic in the central region of the 
trap. We computed the gi(p;0) and g2(p;0) correlations 
functions, which have anisotropic character in the cen- 
ter of the gas even though the density is isotropic in the 
central region of the DBEC. This has measurable implica- 
tions for interference properties and density fluctuations 
of the gas. We found that the anisotropic correlations 
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occur due to roton like quasiparticle modes that can lead 
to local collapse along the polarization axis. 



An immediate extension of this work is to study the 
static structure factor such as was recently done in Ref. 
[22] at zero temperature. An important difference is 
that the static structure will become anisotropic, i.e. 
S(k x , k y ), and with present methods, it could be done at 
finite temperature [41]. This also presents a promising 
method to observe the anisotropic coherences predicted 
here. 



FIG. 2: (a) The Bogoliubov de Gennes spectrum is shown as 
a function of particle number for gd = 0.025 for a/n = 0.25. 
(b,c) A contour plot of the roton quasiparticle modes (u 7 ) 
which go soft for a dipolar system with tilted polarization axis 
are also shown. The dashed black line is »/no an d for w 7 the 



shaded regions are less than zero. The contours at 0.25, 0.5, 
and 0.75 the maximum values of the individual wavefunctions. 




x,y [trap units] 



FIG. 3: (a) The total density (b) the compressibility, kTn(p), 
for the dipolar gas along the x (solid) and y (dashed) axis for 
particles number 2000 (blue), 3000 (red), and 4000 (black). 



Future work will be to investigate the impact of corre- 
lations at higher temperatures and strongly interactions 
where they might impact the BKT transition [U [5] [35] . 
This will require an improvement to the method |40j or 
new approach altogether. Either way, it will be intriguing 
to study the impact of the anisotropic dipolar interaction 
on vortex correlations near the BKT transition and the 
time dependence of such a gas. The present work was 
not able to measure the decay of the correlation func- 
tions due to the relatively weak interaction strength or 
small number the numerical method could handle. If the 
next method can probe stronger interactions, we would 
hope to observe a different character in the decay of the 
correlation functions along and perpendicular to the po- 
larization axis. 
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